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HE: [RH] AMZER E Attagenus unicolor japonicus 的 线 KR EE ZRZEMUA XR RIJ RRRA, 
【方法 ] 本 研究 利用 下 一 代 测 序 方法 测定 了 黑 毛 皮 吉 的 近 完 全 线粒体 基因 组 序列 ,并 基于 39 40 
k H coxl 基因 序列 ,通过 最 大 似 然 法 和 贝 叶 斯 法 构建 了 皮 客 科 的 系统 发 育 树 。 【结果 】 黑 毛皮 过 线 
粒 体 基因 组 包含 37 个 基因 (2 个 rRNA 基因 ,22 个 tRNA 基因 ,13 个 蛋白 质 编码 基因 ) 和 一 段 不 完 
全 的 控制 区 ,总 长 度 为 14 793 bp( GenBank 登录 号 : MG017450) 。 基 因 排 列 顺序 与 已 经 发 表 的 大 多 
数 昆 虫 线粒体 基因 组 相同 ,不 具有 基因 重 排 现 象 。13 个 蛋白 质 编码 基因 中 除 nadl 和 nad2 之 外 其 
他 基因 的 起 始 密码 子 均 是 典型 的 ATN ,而 mnadl 和 nad2 推断 的 起 始 密码 子 是 TTG; cox2, nad5 及 
nad6 基因 终止 密码 子 为 不 完全 的 终止 密码 子 T 外 ,其 他 蛋白 质 编码 基因 的 终止 密码 子 是 典型 的 
TAA 或 TAG。 除 了 trnSl 之 外 ,所 有 的 tRNAs 均 可 形成 典型 的 三 叶 草 结构 ,trnSl 缺少 DHU 辟 而 形 
成 了 一 个 简单 的 环 ,不 能 形成 稳定 的 三 叶 草 结构 。irnSl1 的 反 密 码 子 变 成 了 UCU, 而 不 是 更 常见 的 
GCU, 16S rRNA 的 二 级 结构 包括 5 个 结构 域 (第 I ,全 , MV,，V 和 VW 结构 域 ) ,第 亚 结 构 域 普 遍 缺 
失 而 变 成 一 条 单一 的 核 普 酸 链 连 接 第 卫 和 第 信 结 构 域 , 共 44 个 螺旋 。12S rRNA 的 二 级 结构 包括 4 
个 结构 域 , 共 27 个 螺旋 ,其 中 螺旋 HAT 序列 变异 性 较 高 ,由 一 个 长 共和 一 个 大 环 组 成 。 系 统 发 育 
分 析 结 果 显 示 , 皮 吉 科 内 亚 科 级 别 的 系统 发 育 关系 为 :(( 长 皮 圳 亚 科 (Megatominae) + £  d& 3E 
( Attageninae) ) + X & 3E £f ( Dermestinae) ) ; [2536] X & 3E £F ( Dermestinae) , X & Æ Dermestes 及 
E & k Anthrenus 35] 29 3E R 8E, BE KC M Trogoderma 7; 2 Rž, E € JE ff Dermestinae fe d 44A 
余 的 类 群 构成 姐妹 群 。 
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Analysis of the mitochondrial genome of Attagenus unicolor japonicus 


( Coleoptera: Dermestidae) and a phylogenetic analysis of Dermestidae 
LIN Ai-Li, LI Xin-Xin, ZHAO Xin-Cheng, SONG Nan" ( College of Plant Protection, Henan 
Agricultural University, Zhengzhou 450000 , China) 

Abstract: [ Aim] To analyze the mitochondrial genome of Attagenus unicolor japonicas and the phylogeny 
of Dermestidae. [Methods] The nearly complete mitochondrial genome of A. unicolor japonicus was 
sequenced by next-generation sequencing ( NGS) method, and the phylogenetic trees were constructed 
using maximum likelihood method and Bayesian inference based on coxl gene sequences of 39 insect 
species. [Results] The mitochondrial genome of A. unicolor japonicus contains 37 genes (2 ribosomal 
RNA genes, 22 tRNA genes and 13 protein-coding genes) and a partial control region , with the length of 
14 793 bp (GenBank accession no. : MG017450). No gene rearrangement occurs in the genome 
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organization of this insect, like most Coleoptera mitogenomes published. All protein-coding genes have 
typical start codon ATN except for nadl and mad2 which start with the putative start codon TTG. cox2, 
nad5 and nad6 genes adjacent to the tRNAs terminate with T, while the other protein-coding genes stop 
with the typical stop codon TAA/TAG. The secondary structures of all tRNAs are of typical cloverleaf 
structure, except that of trnS1, in which the dihydrouracil arm ( DHU arm) can not form a stable stem- 
loop structure but forms a simple loop. In addition, the anticodon for trn$1 is UCU rather than the usual 
GCU. The secondary structure of 16S rRNA contains five domains (domain I , IL, IV, V and M) but 
lacks domain lll, with 44 helixes in total. The secondary structure of 12S rRNA consists of four domains 
and 27 helixes. The helix H47 is a highly variable region, which composes of a long stem and a big loop. 
The results of phylogenetic analysis showed that the relationship of Dermestidae subfamilies is 
( ( Megatominae + Attageninae) + Dermestinae). [Conclusion] Dermestinae, Dermestes and Anthrenus 
are all monophyletic, while Trogoderma is polyphyletic. Dermestes is a sister group to the rest of 
Dermestidae. 


Key words: Dermestidae; Attagenus unicolor japonicus; rRNA; secondary structure; mitochondrial 


genome; phylogeny 





BE TEPC Dermestidae ) I 8138] H RE REHPR— 
个 重要 的 具有 检疫 意义 的 储藏 物 害 虫 类 和 群 
( Lawrence and Britton, 1991; Kiselyova and Mchugh, 
2006; Bocak et al., 2014) ,有 6 ^P MEE: K ia MV RE 
( Dermestinae) , [ii] B $8 ME. SF ( Thorictinae ) , J'6 Rz Z 
亚 科 ( Orphilinae ) , E Hz 2€ Mff] ( Trinodinae) ,*E Hz 2k 
亚 科 ( Attageninae ) , K Hz 2 E P} ( Megatominae ) 
( Bouchard et al., 2011) , 4 34 个 属 ,大 约 1 000 种 
( Mroczkowski, 1968; Beal, 1991; Hava, 2003 ) ,我 
国 已 知 有 8 属 40 余 种 。 

皮 埠 科 昆 虫 分 布 较 广 ,食性 复杂 ,主要 危害 动 
物 尸 体 、 各 种 药材 、 毛 织 服装 类 、 皮 张 等 储藏 物 ( 赵 
养 昌 等 , 1966; 张 生 芳 等 , 2007, 2008) ,因此 对 皮 
震 科 昆虫 的 系统 发 育 学 研究 从 经 济 角度 具有 重要 
意义 。 昆 虫 的 线粒体 基因 组 是 一 个 环 状 闭合 的 双 
链 DNA 分 子 ,一 般 包 含 37 个 基因 一 一 13 个 蛋白 
质 编码 基因 (coxl，cox2 ，cox3 cytb, atp6, atp8, 
nadl , nad2, nad3, nad4, nad4L, nad5 和 nad6) 22 
个 转运 RNA (tRNA) 基因 和 2 个 核糖 体 RNA (128 
rRNA 和 16S rRNA) 基因 ,还 有 一 段 控 制 区 (control 
region, CR) ,由 于 AAT 碱 基 含量 较 高 ,也 叫做 A +T 
富 含 区 。 

近年 来 ,对 皮 震 科 的 研究 较 少 , 主要 集中 在 依据 
传统 形态 特征 和 分 子 数据 进行 分 类 鉴定 及 利用 形态 
特征 进行 系统 发 育 研 究 ( Kiselyova and Mchugh, 
2006; KEF S, 2007; BRIF, 2009; 刘 聪 等 ， 
2012; 邓 海 娟 等 , 2014; 周 萍 等 , 2015). HB CL 
止 到 2017 年 9 月 ) GenBank. E E PEZE HS Z E PHI ZR 






































粒 体 基因 组 仅 4 个 ,而 且 均 为 不 完全 的 线粒体 基因 
组 ,已 提交 的 单个 coxl 基因 虽然 有 507 个 ,但 也 是 
不 全 的 基因 片段 且 重 复 较 多 。 本 研究 测定 了 皮 埠 科 
黑 毛 皮 老 Attagenus unicolor japonicus 近 线 粒 体 全 基 
因 组 ,并 对 线粒体 基因 组 进行 了 全 面 分 析 , 基于 
NCBI 上 已 有 的 coxl 基因 3 ' 端 序列 进行 了 皮 霸 科 系 
统 发 育 分 析 , 为 皮 完 科 的 分 类 情况 及 皮 完 科 的 线 粒 
体 基 因 组 研究 提供 科学 数据 。 




















1 材料 与 方法 


1.1 标本 采集 和 DNA 提取 

黑 毛 皮囊 标本 在 2015 年 5 月 采集 于 河南 省 信阳 
鸡 公 山 自 然 保护 区 ,然后 浸泡 在 无 水 乙醇 中 -20Y fg 
存 。 根 据 天 根 公司 的 血液 /细胞 /组 织 基因 组 DNA 提 
取 试 剂 盒 ( TIANamp Genomic DNA Kit) 步骤 提取 黑 
毛皮 圳 昆虫 的 足 和 前 胸 背 板 的 肌肉 获得 总 DNA。 
1.2 线粒体 基因 组 的 获得 

将 提取 的 黑 毛 皮 喜 总 DNA 利用 紫外 分 光 光 度 
YEAR 196 的 琼脂 糖 凝 胶 电 泳 进行 浓度 检测 和 纯度 检 
测 , 检 测 合 格 后 送 与 上 海 欧 易 公 司 利 用 Humina 
HiSeq2500 平台 进行 二 代 测 序 , 共 得 到 约 1 Gb 的 双 
端 测序 长 度 为 150 bp 的 片段 ,将 得 到 的 原始 序列 删 
除 连接 物 .未 配对 的 序列 以 及 一 些 短 的 、 质 量 差 的 序 
列 后 ,高 质量 的 序列 利用 SOAPdenovo-Trans ( Xie et 
al., 2014) 进行 组 装 拼接 。 
1.3 线粒体 基因 组 的 拼接 ,矫正 

我 们 根据 组 装 得 到 的 contigs 利用 Bioedit 7. 0.5.3 
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版 本 (Hall ,1999 ) 进行 本 地 建 库 , 然 后 利用 GenBank 
上 已 提交 的 鞘翅 目 皮 完 科 线粒体 基因 组 的 cox , 
128 rRNA 和 cytb 序列 进行 搜索 比 对 ,接着 把 确认 为 
线粒体 基因 组 的 3 条 序列 上 传 至 Mitos( Bernt et al., 
2013) 进行 注释 分 析 , 将 Mitos 分 析 确 定 的 线粒体 基 
因 组 利用 NCBI Blast ( https: // blast. ncbi. nlm. nih. 
gov/ Blast. cgi? PROGRAM = blastn&PAGE, TYPE = 
BlastSearch&LINK | LOC = blasthome ) 和 BOLD 数据 








JÆ ( http: // www. boldsystems. org/index. php/IDS _ 





OpenldEngine) JETT4) T AAE , 3E 4E 3 RW AE Je: BER 
THE CREE FERE Z EJB Eb do, Aa rub Y RRES 
鉴定 ( 祝 长 清 等 ,1999 ) 。 

根据 GenBank HP E fry Rz ai FIT Zk pr ps di 
组 将 Mitos THE FEE BS CE e ERR XUI DS E DAT 2H. rog 
基因 位 点 进行 手动 矫正 ,其 中 蛋白 质 编码 基因 利用 
MEGA 5.0(Tamura et al., 2011) 进行 翻译 比 对 , 确 
定 基 因 正 确 起 始 、 终 止 位 点 ;tRNAs 和 rRNAs 利用 
MAFFT Version 7 在 线 软件 (https: // mafft. cbrc. jp/ 
alignment/server/ ) 进行 比 对 ,最 终 确定 每 个 基因 的 
正确 起 始 .终止 位 点 。 从 组 装 分 析 结 果 来 看 DUE TN 
区 不 完整 。 
1.4 线粒体 基因 组 分 析 

本 研究 利用 OAGRAW Version 1. 1 (http: // 
ogdraw. mpimp-golm. mpg. de/cgi-bin/ogdraw. pl) 4X 
件 绘制 了 黑 毛皮 震 的 线粒体 基因 组 结构 图 。 

利用 MEGA 5.0(Tamura et al., 2011) 对 黑 毛 皮 
过 线 粒 体 基因 组 的 13 个 蛋白 质 编 码 基 因 进 行 了 碱 
基 组 成 .密码 子 使 用 情况 及 密码 子 使 用 频率 分 析 。 
22 个 tRNA 二 级 结构 依据 Mitos 分 析 结 果 进 行 绘制 。 
黑 毛 皮 圳 的 12S rRNA 和 16S rRNA 的 二 级 结构 依 
据 西方 蜜蜂 Apis mellifera ( Gillespie et al., 2006) fil 
RR Ud Drosophila melanogaster ( Kobayashi and 
Okada, 1990) 的 rRNAs 的 二 级 结构 模型 进行 绘制 。 
1.5 系统 发 育 分 析 
1.5.1 coxl 矩阵 碱 基 蔡 换 饱 和 性 检测 :本 研究 利用 
DAMBE( Xia, 2013) 选择 无 脊椎 动物 线粒体 DNA. 
(invertebrate mitochondrial DNA, InvMtDNA ) 45$ 5 E 
密码 子 表 对 coxl ELCHE BRE PU ETIA H RES AR A 
性 分 析 。 
1.5.2. 系统 发 育 树 构建 :该 研究 基于 已 测 得 黑 毛 皮 
吉 线 粒 体 基因 组 中 coxl 序列 结合 GenBank 中 已 有 
coxl 的 5' 端 序列 (长 度 约 700 bp)39 个 (涉及 34 个 
皮 震 科 物 种 及 作为 外 群 的 5 个 窃 窟 科 物 种 ) 组 成 的 
和 矩阵 利用 最 大 似 然 法 和 贝 叶 斯 法 进行 系统 发 育 构 












































建 。coxl 数据 集 利 用 MEGA 5. 0 ( Tamura et al., 
2011) 进 行 比 对 。 该 研究 利用 IQ-TREE ( Nguyen et 
al., 2015) 网 上 服务 器 (hittp: // iqtree. cibiv. univie. 
ac. at/ ) 基于 coxl 组 成 的 核 昔 酸 序列 矩阵 进行 最 大 
似 然 法 分 析 ,模型 为 CTR ,构建 了 系统 发 育 树 ML- 
cox1DNA ,系统 发 育 树 的 节点 支持 率 由 自 举 检验 置 
信和 度 (3 000 次 重复 ) 进行 评估 。 该 研究 利用 
CIPRES Science Gateway ( http: // www. phylo. org/) 
网 上 服务 器 中 的 Mrbayes 3.2.6 工具 构建 了 贝 叶 斯 
树 BI-cox1 DNA, 其 参数 如 下 :lset nst = 6，rates = 


invgamma，mcmcp ngen = 10 000 000, relburnin = 

















yes, burninfrac = 0. 25, printfreq = 1 000, samplefreq 
-] 000, nruns =2, nchains 2 4, savebrlens = yes, 
节点 支持 率 由 贝 叶 斯 后 验 概 率 ( Bayesian. posterior 
probabilities, PP) 表示。 

















2 结 


2.1 线粒体 基因 组 结构 

黑 毛 皮 震 的 线粒体 基因 组 长 度 为 14 793 pp ,其 
中 A+T 含 量 为 72.38% ,G+C 含量 为 27.62% , 共 
包含 37 个 基因 和 一 段 控 制 区 ,但 控制 区 不 完整 , 仅 
178 bp , 黑 毛 皮 震 线粒体 基因 组 结构 图 如 图 1 所 示 。 
黑 毛 皮 震 的 线粒体 基因 组 注释 如 表 1 Bron ,该 线 粒 
体 基 因 组 已 提交 至 NCBI 数据 库 , 其 登录 号 为 
MG017450, 
2.2. 蛋白 质 编码 基因 分 析 

黑 毛 皮 震 的 线粒体 基因 组 中 13 个 蛋白 质 编码 
基因 总 长 度 为 11 144 bp( 仅 计算 13 个 蛋白 质 编码 
基因 长 度 之 和 ,不 计 基 因 间 隔 ) ,A + 了 含量 为 71%， 
G+C 含 量 为 29% ,其 中 nad3, nad4, nad4L, nad5, 
nad6 ，coxl ，cox2 , cox3 , cytb, atp8 和 atp6 基因 的 起 
台 密 码 子 均 是 典型 的 ATN ,而 nad2 和 nadl 基因 的 
起 始 密码 子 是 TTG; nadl, nad2, nad3, nad4, 
nad4L, coxl, cox3, atp8 和 atp6 基因 的 终止 密码 子 
均 是 典型 的 TAA 和 TAG ,而 cox2, nad5 和 mad6 都 
是 以 T 结尾 。13 个 蛋白 质 编码 基因 序列 的 氨基 酸 
使 用 情况 中 , ESAE Leu) AR (Ser) SEWAR 
(Phe) FEAR (e) 使 用 频率 较 高 ,平均 比例 依次 
为 13.845% , 9. 492% , 9.95196 和 9. 35596 ~ Æ Jt 
AR Cys) FAM (Arg) HAR ( His) RAAR 
(Asp) 使 用 较 低 ,平均 比例 依次 为 1.082% , 1.514% , 
1.839% 和 1. 81296 ~ atp8 中 出 现 的 氨基 酸 种 类 最 
少 , 半 胱 氨 酸 (Cys) RAAR (Asp) 、 谷 氨 酸 (Glu)、 
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TUB BUM 
Attagenus unicolor japonicus 


14 793 bp 








图 1 黑 毛 皮 守 线粒体 基因 组 结构 图 


Fig. 1 Mitochondrial genome structure of Attagenus unicolor japonicus 


加 粗 线 条 外 侧 表 示 基 因 方 向 为 顺 时 针 ,内侧 表 示 基 因 方 向 为 道 时 针 。The outer circle represents genes on the heavy strand (clockwise), while the 





inner circle represents genes on the light strand ( counterclockwise ) . 


甘氨酸 (Gly) 、 组 氨 酸 (His) 和 精 氨 酸 ( Arg ) 都 未 出 
现 ,coxl, nad4L, nadl 和 nad6 序列 中 出 现 氨基 酸 种 
类 约 18 种 ,剩余 蛋白 质 编 码 基 因 序 列 中 出 现 的 氨基 
酸 种 类 均 是 20 种 。13 个 蛋白 质 编码 基因 序列 的 密 
码 子 使 用 情况 如 表 2 所 示 。 
2.3 tRNA 基因 分 析 
REEERE HP 22 个 tRNA 基因 
总 长 度 为 1 435 bp( 仅 计算 22 个 tRNA 基因 长 度 之 
M, AENEA), A +T 含量 为 74.56% ,G+C E 
量 为 25. 44% ,其 二 级 结构 中 除 trasi 之 外 所 有 的 
tRNAs 均 是 典型 的 三 叶 草 结构 ,trnsl 缺少 DHU 臂 











而 形成 了 一 个 简单 的 环 ,不 能 形成 稳定 的 三 叶 草 结 
FJ. 535b .trnS1 的 反 密码 子 变 成 了 UCU ,而 不 是 更 
常见 的 GCU。 碱 基 配 对 中 除了 上 典型 的 Waston-Crick 
配对 A-U.G-C 外 ,G .0U 配 对 也 是 存在 的 , 称 之 
为 非典 型 配对 , 也 称 为 摆动 配对 (wobble base 
pairs) 。 在 22 个 tRNAs 的 二 级 结构 中 , 出现 一 些 碱 
基 错 配 及 非 规范 配对 现象 。:rzW 中 出 现 了 AC 碱 基 
对 错 配 情况 ,trmW 出 现 AG 碱 茜 对 错 配 情况 ,trnLl， 
irnL2 , trnT , trnF 和 trnV 出 现 UU 碱 基 对 错 配 情况 ， 
trnC, trnY, trnG, trnA, trnSl, trnH, trnP, trnF 及 
trnLl 出 现 非 规范 性 碱 基 对 G + U, 
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R1 ， 黑 毛皮 圳 线粒体 基因 组 注释 


Table 1 Annotation of the mitochondrial genome of Attagenus unicolor japonicus 























基因 ”基因 长 度 (bp) 起 始 位 置 (bp) 终止 位 置 (bp) ”基因 间隔 (bp) ”起 始 密 码 子 ”终止 密码 子 MRF 编码 链 
Gene Gene length Start. position Stop position ~ Intergenic length Start codon Stop codon Anticodon Coding strand 
trn] 65 65 129 -3 GAT H 
trnQ 69 127 195 -1 TTG L 
irnM 69 195 263 0 CAT H 
nad2 1 005 264 1 268 2 TTG TAA H 
trnW 65 1271 1 335 -8 TCA H 
trnC 63 1 328 1 390 0 GCA L 
trnY 64 1 391 1 454 -8 GTA L 
coxl 1 548 1 447 2 994 -5 ATT TAA H 
trnL2 66 2 990 3 055 0 TAA H 
cox2 682 3 056 3 737 0 ATA T H 
trnK. 69 3 738 3 806 0 CTT H 
trnD 66 3 807 3 872 0 GTC H 
atp 156 3 873 4 028 -7 ATC TAA H 
atp6 669 4 022 4 690 -1 ATG TAA H 
cox3 789 4 690 5 478 -1 ATG TAA H 
trnG 66 5 478 5 543 0 TCC H 
nad3 351 5 544 5 894 -2 ATA TAG H 
trnA 65 5 893 5 957 0 TGC H 
traR 64 5 958 6 021 3 TCG H 
trnN 65 6 025 6 089 0 GTT H 
trnS1 67 6 090 6 156 0 TCT H 
trnE 63 6 157 6 219 -2 TTC H 
trnF 63 6218 6 280 0 GAA L 
nad5 1714 6 281 7 994 0 ATT T L 
trnH 62 7 995 8 056 -1 GTG L 
nad4 1 332 8 056 9 387 -7 ATG TAA L 
nad4L 282 9 381 9 662 10 ATA TAA L 
trnT 62 9 673 9 734 0 TGT H 
trnP 64 9 735 9 798 1 TGG L 
nad6 525 9 800 10 324 -5 ATT TAA H 
cytb 1 128 10 320 11 447 -2 ATT TAA H 
trnS2 67 11 446 11 512 18 TGA H 
nadl 951 11 531 12 481 1 TTG TAA L 
trnLl 61 12 483 12 543 0 TAG L 
16S rRNA 1 257 12 544 13 800 0 L 
traV 70 13 801 13 870 0 TAC L 
128 rRNA 745 13 871 14 615 0 L 
Partial CR 178 14 616 14 793 64 非 编码 序列 


Non-coding sequence 


H: 重 链 Heavy strand; L: 轻 链 Light strand. T 表示 不 完全 终止 密码 子 ,通常 认为 是 转录 后 的 mRNA 的 3 "端的 A 补充 为 完整 的 终止 密码 子 。“T” 
indicates the incomplete stop codon. It is generally assumed that TAA stop codon is completed by the addition of 3' A to the mRNA. 

















482 昆虫 学 报 Acta Entomologica. Sinica 61 卷 
表 2， 黑 毛皮 圳 线粒体 基因 组 13 个 蛋白 质 编码 基因 的 相对 密码 子 使 用 频率 
Table 2 Relative synonymous codon usage of 13 protein-coding genes of the mitochondrial genome 
of Attagenus unicolor japonicus 
密码 子 。 ”数量 密码 子 。 ”数量 密码 子 。 ”数量 密码 子 。 ”数量 
Codon Count ius Codon Count dna Codon Count ER Codon Count ESQ 
UUU(F) 314 1.48 UCU(S) 45 1.09 UAU(Y) 181 1.57 UGU(C) 43 1.32 
UUC(F) 111 0.52 UCC(S) 33 0. 80 UAC(Y) 50 0.43 UGC(C) 22 0. 68 
UUA(L) 154 1.96 UCA(S) 86 2.08 UAA(*) 135 1.32 UGA(W) 63 1.17 
UUG(L) 83 1.06 UCG(S) 16 0.39 UAG(*) 69 0.68 UGG(W) 45 0.83 
CUU(L) 93 1.18 CCU( P) 42 1.54 CAU(H) 60 1.45 CGU(R) 4 0.38 
CUC(L) 32 0.41 CCC( P) 26 0.95 CAC(H) 23 0.55 CGC(R) 5 0.48 
CUA(L) TI 0.98 CCA(P) 33 1.21 CAA(Q) 73 1.51 CGA(R) 22 2.10 
CUG(L) 33 0.42 CCG( P) 8 0.29 CAG(Q) 24 0.49 CGG(R) 11 1.05 
AUU(I) 243 1.57 ACU(T) 51 1.23 AAU(N) 152 1.46 AGU(S) 39 0.94 
AUC(I) 66 0. 43 ACC(T) 31 0.75 AAC(N) 56 0.54 AGC(S) 29 0. 70 
AUA(M) 134 1.48 ACA(T) 73 1.76 AAA(K) 108 1.58 AGA(S) 42 1.02 
AUG(M) 47 0.52 ACG(T) 11 0.27 AAG(K) 29 0.42 AGG(S) 41 0.99 
GUU( V) 64 1.63 GCU( A) 19 1.03 GAU(D) 58 1.51 GGU(G) 36 0.94 
GUC( V) 21 0.54 GCC(A) 25 1.35 GAC( D) 19 0.49 GGC(G) 13 0.34 
GUA( V) 55 1.40 GCA( A) 26 1.41 GAA(E) 61 1.45 GGA(G) 58 1.51 
GUG( V) 17 0. 43 GCG( A) 4 0.22 GAG(E) 23 0.55 GGG(G) 47 1.22 











RSCU: 相对 密码 子 使 用 频率 Relative synonymous codon usage. 





2.4 rRNA 基因 分 析 

本 研究 依据 西方 蜜蜂 A. mellifera ( Gillespie et 
al., 2006) ) E JR Rta D. melanogaster ( Kobayashi 
and Okada, 1990) 的 rRNAs 的 二 级 结构 模型 绘制 了 
REB EB 12S rRNA 和 16S rRNA 的 二 级 结构 。 
128 rRNA 长 度 为 745 bp ,位 于 troV 和 控制 区 之 间 ， 
A €T ig d 75.3096 ,G +C 含量 为 24.70% ,其 二 
级 结构 包括 4 个 结构 域 ( 第 工 - 玉 结构 域 ) , 共 27 
个 螺旋 ,其 中 螺旋 H47 序列 变异 性 较 高 ,由 一 个 长 
茎 和 一 个 大 环 组 成 。16S rRNA 长 度 为 1 257 pp ,位 



































* 终止 密码 子 Stop codon. 


于 trnLl 和 trnV 之 间 , A + 了 TT 含量 为 78.28% ,G+C 
含量 为 21. 72% ,其 二 级 结构 包括 5 个 结构 域 (第 
L, I, W, V 和 VY 结构 域 ) ,第 亚 结构 域 普 遍 缺 失 
而 变 成 一 条 单一 的 核 苷 酸 链 连接 第 工 和 第 IV 结构 
域 , 共 44 个 螺旋 。 
2.5 coxl 数据 集 的 系统 发 育 信 号 检测 

本 研究 利用 DAMBE( Xia, 2013) 软件 对 coxl 核 
苷 酸 序列 进行 核 苷 酸 蔡 代 饱 和 性 分 析 。 结 果 如 表 3 
所 示 ,39 个 类 群 的 coxl 序列 密码 子 第 3 位 点 出 现 饱 
和 ,第 1 位 点 和 第 2 位 点 比较 保守 。 























表 3 PA DAMBE 检测 coxl 序列 碱 基 蔡 换 饱 和 性 
Table 3 Sequence saturation test of coxl performed by DAMBE 














数据 集 Data base NumOTU Iss Iss. cSym P Iss. cAsym P 
39taxa-coxl 序列 的 第 1 位 点 39taxa-coxl_1 32 0. 184 0. 686 0. 0000 0. 367 0. 0000 
39taxa-coxl 序列 的 第 2 位 点 39taxa-coxl_2 32 0. 095 0. 686 0. 0000 0. 367 0. 0000 
30taxa-coxl 序列 的 第 3 位 点 39taxa-coxl_3 32 0.744 0. 686 0. 0578 0. 367 0. 0000 
39taxa-coxl 序列 的 第 1 和 2 位 点 39taxa-coxl_1 42 32 0. 134 0.695 0.0000 0. 367 0. 0000 
391axa-coxl 的 核 昔 酸 全 序列 39taxa-coxl 1 42 +3 32 0. 290 0.718 0. 0000 0. 392 0. 0000 
NumOTU ; 可 操作 分 类 单元 数量 Number of operational taxonomic units; P; 概率 Probability; Iss: 替换 饱和 指数 Index of substitution saturation; /ss. 
cSym ; 假设 拓扑 结构 对 称 Iss. c assuming a symmetrical topology; Iss. cAsym: 假设 拓扑 结构 不 对 称 Iss. c assuming an asymmetrical topology. Iss < Iss. 


c: 不 饱和 None saturation; Iss > 1ss. c; 替换 饱和 Substitution saturation. 出 现 饱和 的 数据 加 粗 显 示 。Saturation data are indicated in bold. 
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2.6 系统 发 育 分 析 
DERHEN PME 





JET BORERE34 个 种 的 coxl 











d PL SHEER UA ER Fg ADAE RI 11 mp9 1 
构建 的 系统 发 育 树 ML-coxI DNA 和 BI-cox1 DNA 如 


图 2 和 3 所 示 ,两 个 系 乡 





发 育 树 的 拓扑 结构 大 致 相 


同 。 皮 震 科 内 亚 科 级 别 的 系 乡 


发 育 关 系 为 ,(( 长 皮 


34 FE Megatominae) + E ZPH Attageninae ) ) 
+ Æ a MERE ( Dermestinae )), K Æ Z& NV. $ 


( Attageninae ) R Jy — 5 , E Hz 3$ EP} ( Dermestinae ) 
构成 姐妹 群 , C AA SC UE 100, K BE 3X EE 
( Megatominae) 15 SHP EZ R J Wi AKIE : [p Hz dk 
属 Anthrenus 为 单 系 群 ,与 Ctesias PIEK ff [9] Sz a JB 
Globicornis & Jj — fR ; XXE Rz 48 Ja Trogoderma 为 多 系 
群 ,与 Anthrenocerus 和 Reesa ER 3 — Iko "E HERE WV. 
F} ( Attageninae ) 9 个 物种 聚 为 一 文 , E RES 8 








Attagenus 为 并 系 群 。 





F8 MV RE ( Dermestinae ) 10 个 


( Megatominae ) ^j E RE d EPH IRAR RE , 3 e LRE 
率 为 96。 长 皮 喜 亚 科 (Megatominae ) ME Rz 5 MV f 


WIRA — 3c , BERE JR. Dermestes 为 单 系 群 。 








Anthrenus festivus (KU. 494037). 一 
JUCERTRL Re Anthrenus scrophulariae (KM. 442953) 
B S [E] Hz afe Anthrenus flavipes (KU. 494040) 
EE E He d Anthrenus pimpinellae (KM. 441362) 
Tas BL AER. Anthrenus museorum (KM. 447423) 
VQ RES Anthrenus sp. (KR. 130254) 
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BERE SIE Trogoderma glabrum (KM 441814). 一 
EFE Trogoderma anthrenoides (KP. 331485) 
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Sefrania bleusei (KU 494237) 


Attagenus trifasciatus (KU. 494075) — 


伍 氏 毛皮 Attagenus woodroffei (KJ. 965226) 


横 带 毛皮 咱 Attagenus fasciatus (KU 494069) 


WEEE Attagenus smirnovi (KU. 494070) 


15.085 ik Attagenus unicolor (KM_446244) 
Eik Attagenus unicolor japonicus (MG 01745 
暗 褐 毛皮 党 Attagenus brunneus (KU 494058) 
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ZEER Attagenus pellio (KM. 443611) 


WAREZ Dermestes frischii (KP_400413) — 


FTBEHZ 3 Dermestes maculatus (HM. 909035) 
HARE Dermestes undulatus (KU. 494100) 
Dermestes szekessyi (KJ 961991) 
SA HE Dermestes laniarius (KM. 442171) 
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PET Dermestes bicolor (KM. 441916) 
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FRP Z Anobium punctatum (KM. 447362) 


MAT Lasioderma serricorne (KM. 216269) 


Gastrallus immarginatus (KM 443499) 
药材 甲 Stegobium paniceum (KM. 452133). — 


DERNE) 


Anobiidae (Out-group) 


图 2 最 大 似 然 法 构建 的 基于 39 个 物种 的 coxl 核 背 酸 序列 的 系统 发 育 树 
Fig. 2 Phylogenetic tree of 39 species based on nucleotide sequence of coxl using maximum likelihood method 
节点 旁 数字 代表 ML 树 的 支持 率 。Node numbers refer to bootstrap support values. 


3 结论 与 讨论 


3.1 线粒体 基因 组 结构 


黑 毛 皮 寺 的 线粒体 基因 组 长 度 为 14 793 bp, £i 


含 37 个 基因 和 一 段 控制 区 ,其 线粒体 基因 组 的 方向 
及 基因 排列 顺序 都 和 假定 同一 祖先 的 不 同 昆虫 的 线 
粒 体 基 因 组 完全 一 样 (Boore，1999 ) , 2E BERE 
线粒体 基因 组 13 个 蛋白 质 编码 基因 中 , BR cox2 ， 
nad5 及 nad6 终止 密码 子 以 了 结尾 之 外 ， 所 有 和 蛋白 
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斯 法 构建 的 基于 39 个 物种 的 coxl 核 苷 酸 序列 的 系统 发 育 树 


Fig. 3 Phylogenetic tree of 39 species based on nucleotide sequence of coxl using Bayesian inference 


节点 旁 数 字 代 表 BI 树 的 后 验 概率 。Node numbers refer to Bayesian posterior probabilities. 





质 编码 基因 的 终止 密码 子 均 是 典型 的 TAA 和 TAG, 
通常 认为 是 转录 后 的 mRNA 的 399 B A 补充 为 完 
整 的 终止 密码 子 。 

黑 毛 皮 守 的 线粒体 基因 组 22 个 tRNA 基因 中 
TRI irnSl 之 外 所 均 是 典型 的 三 叶 草 结构 ,irnSl dA 
少 DHU 臂 而 形成 了 一 个 简单 的 环 ,不 能 形成 稳定 的 
三 叶 草 结构 ,这 在 昆虫 纲 其 他 类 群 的 昆虫 中 也 出 现 
过 (Bae et al., 2004; Kim et al., 2005; Stewart and 
Beckenbach , 2005 ; Li et al., 2012; Sun et al., 2012; 
Wang et al., 2012; Wang and Tang, 2017) 。 另 外 ， 
trnSl 的 反 密 码 子 变 成 了 UCU 而 不 是 更 常见 的 
GCU, 这 种 情况 在 畏 翅 目 甲 虫 研究 中 也 有 出 现 
( Stewart and Beckenbach, 2003; Friedrich and 
Muquim, 2003; Sheffield et al., 2008; Wang et al., 











2012) , Friedrich 和 Muquim ( 2003 ) 认为 trnSl 的 反 
密码 子 由 GCU 变 为 UCU 在 节肢 动物 的 演化 中 多 次 
出 现 的 (Hickerson and Cunningham，2000; Shao et 
al., 2001) Sheffield 等 (2008 ) 研究 表明 ,虽然 节肢 
动物 中 trnSl 的 反 密 码 子 一 般 是 GCU ,但 保守 的 反 
密码 子 UCU 在 多 食 亚 目 中 可 能 是 一 个 共 源 性 状 。 
HEH 12S rRNA 二 级 结构 中 ,螺旋 HAT 
序列 变异 性 较 高 ,不 同 模型 的 二 级 结构 也 不 相同 , 黑 
TREES RE HAT 结构 由 一 个 长 蕉 和 一 个 大 环 组 
成 ,这 和 Gillespie 等 (2006 ) 绘制 的 西方 蜜蜂 A. 
mellifera 的 12S rRNA 二 级 结构 中 的 螺旋 H47 结构 
一 致 ,但 Zhang 等 (2016 ) 描述 的 黄粉 虫 Tenebrio 
molitor 的 12S rRNA 二 级 结构 中 的 螺旋 H47 则 是 由 
2 个 长 茎 .1 个 短 茎 和 3 个 环 组 成 。 在 16S rRNA 二 
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级 结构 中 包括 5 个 结构 域 (第 1, 卫 , V,，V 和 VY 结 
构 域 ) ,第 焉 结构 域 通常 缺失 变 成 一 条 单一 的 核 苷 
酸 链 将 第 开 结 构 域 和 第 及 结构 域 连 接 到 一 起 。 
3.2 系统 发 育 分 析 

根据 最 大 似 然 法 获得 的 ML-coxI DNA 树 和 贝 叶 
斯 法 获得 的 BI-cox DNA 树 拓扑 结构 一 致 ,系统 发 
ADRK EEP ERREF KARHE 
视 为 一 个 单 系 群 (Lawrence and Newton, 1982; 
Beutel, 1996; Zhantiev, 2000) , Kiselyova 和 Mchugh 
(2006) XT 4J] EE SRE EE Be AUREAS Dt T x 
8j, IPEP ER AR Orphilus 后 是 单 系 群 ,但 本 研究 
cox 数据 集中 并 未 有 Orphilus 物种 的 数据 参与 。 皮 
A8 FE VJ MERE RU B3 REES OR IR LK BERE ERI 
( Megatominae ) 5j Æ Hz Z MV. E. ( Attageninae ) 构成 姐 
TRE, K Hz 3$ E P ( Megatominae ) F0 Æ Hz Za E P 
( Attageninae ) R Jy — 3c Ej Rz 38 WP} ( Dermestinae ) 构 
成 姐妹 群 。Zhantiev(2000 ) 3E F Bz $8 PHR R 1 A 
REN BRES RIE E T BERRREAS EAE ,该 
WMR A A mA A E BRE NER 
( Dermestinae ) ,J'5 Hz ££ MV EE ( Orphilinae ) X [Ei] Hz za V. 
科 ( Anthreninae ) ,Zhantiev(2000) 系 统 进 化 树 结果 与 
本 文系 统 发 育 结果 几乎 一 致 ,都 支持 皮 寇 亚 科 
( Dermestinae ) 2j !& Z& At, H xe d Hz xà WE 
( Dermestinae) +5 KE 45 PIRI Zi Z8 RE F9 R AH b HE o 
Kiselyova 和 Mchugh( 2006 ) 1H FX 2i Rth ze Hj] Hz 38 MV. 
科 (Dermestinae ) 5 RE SE Fl 2S EF AH IKE, 1E 
是 由 于 Zhantiev (2000 ) 1 X 3 BEA) 2g 3 个 亚 科 , E 
B $E M Ph ( Attageninae ) Wü K Ez 2X XE ORB 
( Megatominae) 在 其 系统 发 育 结 果 中 未 有 命名 ,但 是 
从 属 的 级 别 上 来 看 , 现在 被 划分 为 毛皮 寇 亚 科 
( Attageninae ) FI K Hz 2£ MF £F ( Megatominae ) 下 的 属 
也 聚 为 一 支 , 这 点 和 本 研究 的 系统 发 育 树 结果 是 一 
致 的 ,但 由 于 本 研究 取样 仅 包括 3 P IURE, EE BERE 
科 (Megatominae ) 和 毛皮 圳 亚 科 (Attageninae ) 构成 
姐妹 群 的 关系 还 有 待 其 他 证 据 的 补充 和 验证 。 本 研 
究 系统 发 育 树 的 总 体 节 点 支持 率 较 高 ,但 个 别 节 支 
位 点 支持 率 较 低 ,原因 可 能 是 取样 较 少 ,有 些 属 未 取 
到 样 ,有 些 属 仪 取 一 个 种 ,这 是 由 于 皮 堵 科 物 种 中 已 
有 的 coxl 序列 还 不 够 全 面 ; 男 外 ,一 般 情况 下 , 越 靠 
近 树 的 根部 节点 支持 率 会 降低 ,但 本 研究 构建 的 系 
统 发 育 树 的 根部 节点 支持 率 较 高 ,分 支 底部 反而 降 
低 , 这 是 因为 本 研究 选取 的 皮 老 科 的 coxl 5' 端 序列 
部 分 物种 序列 一 致 性 较 高 ,但 并 不 是 序列 一 致 性 越 
高 ,其 节点 支持 率 就 越 高 。 当 序列 一 致 性 过 高 时 ,其 




































































































































































序列 差异 度 更 小 ,不 易 区 分 ,这 类 序列 与 其 他 序列 进 
行 构 建 系统 发 育 树 的 时 候 聚 在 一 起 的 可 能 性 会 相 
同 , 其 构建 的 系统 发 育 树 就 会 出 现 偏差 ,节点 支持 率 
也 会 降低 。 皮 窟 亚 科 (Dermestinae ) VJ B JZ 28 JE 
Dermestes KK K aa EP} ( Megatominae ) Pj 83 [8] Sz d 
属 Anthrenus JER Z& HE , E EE Jg Trogoderma 是 多 
系 群 这 一 系列 结论 与 Kiselyova 和 Mchugh ( 2006 ) jf 
究 结果 一 致 。 

Castalanelli 等 ( 2011 ) 利用 线粒体 基因 cox , 
cytb 和 18S rRNA FOEDE HERERHIER RART A 
Vr, RE HH DE HZ 33S Trogoderma 是 多 系 群 ,并 估 
43 PREMAZA 1. 75 亿 年 。 邓 海 娟 等 
(2014) EF Ez SEE 3 属 的 线粒体 基因 组 的 coxl FE 
列 进行 系统 发 育 分 析 表 明 ,线粒体 基因 组 中 的 coxl 
基因 可 以 作为 分 子 水 平 上 的 依据 进行 皮 赛 科 昆 虫 生 
物 系统 发 育 分 析 。 本 研究 虽然 解释 了 皮 宕 科 中 一 部 
分 类 群 的 系统 发 育 关系 ,但 取样 仍 较 少 , 且 仅 有 部 分 
coxl 序列 ,分子 依 据 较 少 。 因 此 ,未 来 的 研究 应 该 加 
大 皮 完 科 的 取样 物种 数量 和 线粒体 基因 组 的 数量 ， 
利用 更 多 数据 去 分 析 皮 埠 科 的 系统 发 育 及 皮 才 科 的 
分 类 地 位 。 
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